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Abstract 

The spin network simulator model represents a bridge between 
(generalized) circuit schemes for standard quantum computation and 
approaches based on notions from Topological Quantum Field Theo- 
ries (TQFT). More precisely, when working with purely discrete uni- 
tary gates, the simulator is naturally modelled as families of quan- 
tum automata which in turn represent discrete versions of topolog- 
ical quantum computation models. Such a quantum combinatorial 
scheme, which essentially encodes SU (2) Racah-Wigner algebra and 
its braided counterpart, is particularly suitable to address problems 
in topology and group theory and we discuss here a finite states- 
quantum automaton able to accept the language of braid group in 
view of applications to the problem of estimating link polynomials in 
Chern-Simons field theory. 
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1 Introduction 



In the past few years there has been a tumultuous activity aimed at in- 
troducing novel conceptual schemes for quantum computing. The approach 
proposed in P, [2j relies on the (re)coupling theory of SU{2) angular mo- 
menta (Racah-Wigner tensor category) and might be viewed as a general- 
ization to arbitrary values of the spin variables of the usual quantum circuit 
model based on qubits and Boolean gates. Computational states belong to 
finite-dimensional Hilbert spaces labelled by both discrete and continuous 
parameters, and unitary gates may depend on quantum numbers ranging 
over finite sets of values as well as continuous (angular) variables. When 
working with purely discrete unitary gates, the computational space of the 
simulator is naturally modelled as families of quantum automata which in 
turn represent discrete versions of topological models of quantum computa- 
tion based on modular functors of SU{2) Chern-Simons theory (c/r. jSj and 
references therein). The discretized quantum theory underlying our scheme 
actually belongs to the class of SU{2) 'state sum models' introduced in pi 
and widely used in 3-dimensional simplicial quantum gravity (c/r. section 5 
of [2j and references therein). From the computational viewpoint, we are in 
the presence of finite-states and discrete-time machines able to accept any 
(quantum) language compatible with the Racah-Wigner algebra on the one 
hand, and as powerful as Freedman's quantum field computer on the other. 
As is getting more and more evident, the exponential efficiency that quantum 
algorithms may achieve vs. classical ones might prove especially relevant in 
addressing problems in which the space of solutions is not only endowed with 
a numerical ('digital') representation but is itself characterized by some ad- 
ditional 'combinatorial' structure, definable in terms of the grammar and the 
syntax of some 'language'. Our scheme can be easily adapted -by braiding 
the Racah-Wigner tensor category as usually done in SU{2) Chern-Simons 
framework- to make it accepting the language of the braid group -B„ and 
to handle with (colored) link polynomials expressed as expectation values of 
composite Wilson loop operators. In the last section we are going to exhibit 
a quantum automaton calculation which processes efficiently -linearly in the 
number of pairwise braidings- the language of the braid group. This does not 
mean that we have really got a 'quantum algorithm' since we should show 
that each 'elementary' (in the sense of the spin network simulator) unitary 
gate can be evaluated (or approximated) efficiently. This could be achieved 
by exploiting both the combinatorial properties of the spin network graph 
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and suitable recursion relations which hold for hypergeometric-type polyno- 
mials of discrete variables (work is in progress in such directions). 
Finally, let us point out that the idea of using braiding operators to imple- 
ment quantum gates actually dates back to Freedman and collaborators jS] 
and has been exploited also by Kauffman and Lomonaco [S]. The common 
challenge of all approaches is, on the one hand, the search for suitable unitary 
representations of the braid group and, on the other, the selection of suit- 
able encoding maps into quantum circuits (or, eventually, into new quantum 
computing schemes). Aharonov, Jones and Landau have recently provided in 
jH] an efficient quantum algorithm which approximates the Jones polynomial 
built from Markov traces arising from representations of the braid group in 
the Temperly-Lieb algebra. 

2 Complexity of braids and links 

Let us recall some basic properties of the Artin braid group Bn has n 
generators, denoted by {ai,a2, ■ ■ ■ ,cr„_i} plus the identity e, which satisfy 
the relations 

aiCTj = aj ai if |i — j | > 1 {i, j = 1,2, . . . ,n — 1) 

CTi cn+i = cTi+i CTi+i ( i = 1, 2, . . . , n - 2). (1) 

An element of the braid group is a 'word' w in the standard generators of 
Bn, e.g. w = a^^a2 cr^^a2crf cr2"^(Ticr^^ G B^; the length \w\ of the word w is 
the number of its 'letters'. The group acts naturally on topological sets of n 
disjoint strands - running downward and labeled from left to right - in the 
sense that each generator corresponds to the over-crossing of the ith strand 
on the (z + l)th, and cr,j"^ represents the inverse operation (under-crossing) 
according to cTj a^^ = u^^ai = e. 

As is well known, representation theory of Artin braid group enters heavily 
into many physical applications, ranging from statistical mechanics to (topo- 
logical) quantum field theories. Motivated by this remark, the search for 
algorithms addressing computational problems arising in braid group con- 
text is becoming more and more compelling. Historically, three fundamental 
decision problems for any finitely presented group G were formulated by Max 
Dehn in 1911: 

• word problem: does there exist an algorithm to determine, for any arbitrary 
word w in the generators of G, whether or not w = identity in G? 
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• coniugacy problem: does there exist an algorithm to decide whether any 
pair of words in the generators of G are conjugate to each other? 

• isomorphism problem: given an arbitrary pair of finite presentations in 
some set of generators, does there exist an algorithm to decide whether the 
groups they present are isomorphic? 

Following the development of the classical theory of algorithms (recur- 
sive functions and Turing machine) it is reasonable to expect that Dehn's 
problems might be recursively solvable or, at least, that some 'local' ones 
(the word and the coniugacy problems) be so. It turns out, instead, that 
not only these problems, but a host of local and global decision problems 
sharing a combinatorial flavor are unsolvable within such scheme. As for the 
braid group, and referring to |7j (which contains an exhaustive review to- 
gether with an up-to-date bibliography) , the word problem has been recently 
solved. More precisely, given two words w,w' G Bn there exists an efficient 
algorithm whose time complexity function (namely the time required to per- 
form the computation as a function of the input 'size') grows as Iwl"^ nlog n, 
where \w\ is the length of the word w and n the index of the braid group. 
Thus this problem belongs to the (classical) complexity class P, which con- 
tains languages accepted by a Turing machine in polynomial time. The word 
problem seems to be one short step away from the non-minimal braid prob- 
lem: given a word w in the generators and their inverses, determine 
whether there is a shorter word w' in the same generators which represents 
the same element of Bn- The decision process can be described as follows: 
list all the words that are shorter than the given one, drop from the list as 
many candidates as possible by simple criteria, and then test, one by one, 
whether the survivors represent the same element as w. Surprisingly enough, 
this problem turns out to be NP-complete, namely belongs to the complex- 
ity class of Non deterministic Polynomial algorithm, for which a proposed 
solution can be checked efficiently (polynomially). NP-complete problems 
are in NP and can be characterized by saying that a polynomial algorithm 
for one of them provides automatically an efficient solution for all NP prob- 
lems. Finally, the best known classical algorithm for the coniugacy problem 
in the braid group is exponential in both n and \w\. 

Knot theory is closely related to braid group owing to Alexander's The- 
orem, which states that every knot (or link ) L in the 3-sphere can be 
represented (not uniquely) as a closed braid for some suitable n. Moreover, 
the reduction from a planar diagram of the link L to the closed braid i&(„)(L), 
with W(^n) £ Bn, has been recently shown to be polynomial in the number of 
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strands (see the original references in 0). The major challenge in knot the- 
ory is of course the detection and classification of all possible knots and links 
up to ambient isotopy. There are a number of topological, combinatorial and 
algebraic invariants, starting from the group of the knot (the fundamental 
group of the complement of the knot in S^) and numerical invariants (linking 
number, crossing number, ..) up to invariants of polynomial type. This last 
class of invariants is related to the braid group, in the sense that knot poly- 
nomials are actually associated with characters (technically, Markov traces) 
of (suitable) representations of Bn for some n. A crucial role is played by 
the Jones polynomial [H] , which, when evaluated at particular roots of unity, 
turns out to be associated with suitable expectation values of observables in 
Chern-Simons topological field theory [0]. From the computational side, it 
has been proved that the exact evaluation of the polynomial Vl{uj) at a; = 
root of unity can be performed in polynomial time in terms of the number 
of crossings of planar diagram of L if a; is a 2nd, 3rd, 4th, 6th root of unity. 
Otherwise, the problem is #P-hard [TU] (the computational complexity class 
T^P-hard is the enumerative analog of NP-complete problems). 

Generally speaking, the exponential efficiency that quantum algorithms 
may achieve vs. classical ones might prove especially relevant in addressing 
problems in which the space of solutions is not only endowed with a numerical 
representation but is itself characterized by some additional 'combinatorial' 
structure, definable in terms of a grammar and a syntax and thus suitable to 
be encoded naturally in the spin network computing framework (which will 
be addressed in the following). There are a number of problems that are not 
easily formulated in numerical terms and that are quite often intractable in 
classical complexity theory (c/r. PI]). In combinatorial and algebraic topol- 
ogy typical issues are: the construction of presentations of the fundamental 
group (or the first homology group) of compact 3-manifolds decomposed as 
handlebodies; the study of equivalence classes of knots/links in the three- 
sphere, related in turn to the classification of hyperbolic 3-manifolds; the 
enumeration of inequivalent triangulations of D-dimensional compact man- 
ifolds. 

Of course efficient solutions of this kind of problems would be interesting also 
in view of applications to (discretized) models for quantum gravity at least 
in low spacetime dimensions. 
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3 Spin network quantum simulator 



The theory of binary couphng of = n + 1 SU (2) angular momenta rep- 
resents the generahzation to an arbitrary of the couphng of two angular 
momentum operators Ji, J2 which involves Clebsch-Gordan (or Wigner) co- 
efficients in their role of unitary transformations between uncoupled and 
coupled basis vectors, \jimi > ®|j2'^2 > and \jij2]JM > respectively. 
The quantum numbers ji,j2 associated with Ji, J2 label irreducible repre- 
sentations of SU{2) ranging over {0, 1/2, 1, 3/2, . . .}; mi, 777,2 are the magnetic 
quantum numbers, —ji < rrii < ji in integer steps; J is the spin quantum 
number of the total angular momentum operator J = Ji + J2 whose mag- 
netic quantum number is M = 777i -|- 777.2, — J < M < J. Here units are 
chosen for which h = 1 and we refer to ^2] for a complete account on the 
theory of angular momentum in quantum physics. On the other hand, SU{2) 
'recoupling' theory -which deals with relationships between distinct binary 
coupling schemes of angular momentum operators- is a generalization to 
any of the simplest case of three operators Ji, J2, J3 which calls into play 
unitary transformations known as Racah coefficients or 6j symbols. A full 
fledged review on this advanced topic in the general framework of Racah- 
Wigner algebra can be found in [T^ . 

The architecture of the 'spin network' simulator worked out in (2] relies 
extensively on recoupling theory and can be better summarized by resorting 
to a combinatorial setting where the computational space is modelled as an 
5'f/(2)-flber space structure over a discrete base space V 

{V, C''+\ SU{2)\ (2) 

which encodes all possible computational Hilbert spaces as well as unitary 
gates for any flxed number A^ = 77 + 1 of incoming angular momenta. 

• The base space V = {f (b)} represents the vertex set of a regular, 
3-valent graph Gn{V,E) whose cardinality is \V\ = {2n)\/n\. There exists a 
one-to-one correspondence 

{v{h)}^{ni{h)} (3) 

between the vertices of Gn{V,E) and the computational Hilbert spaces of 
the simulator. 

The label b above has the following meaning -on which we shall exten- 
sively return later on: for any given pair (77, J), all binary coupling schemes 
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of the n + 1 angular momenta {J^}, identified by the quantum numbers 
ji, . . . ,jn+i plus ki, . . . , kn-1 (corresponding to the n — 1 intermediate an- 
gular momenta {Kj}) and by the brackets defining the binary couplings, 
provide the 'alphabet' in which quantum information is encoded (the rules 
and constraints of bracketing are instead part of the 'syntax' of the resulting 
coding language). The Hilbert spaces TC^ {ki, . . . , thus generated, each 
(2 J + l)-dimensional, are spanned by complete orthonormal sets of states 
with quantum number label set B such as, e.g. for n = 3, {((ji (72^3)^1 )fc2^4)j 

More precisely, for a given value of n, 7Y;((b) is the simultaneous eigenspace 
of the squares of 2[n+l) Hermitean, mutually commuting angular momentum 
operators Ji, J2, J3,...,J„+i with fixed sum Ji + J2 + J3 +.. .+Jn+i = J, 
of the intermediate angular momentum operators Ki, K2, K3, . . . , K„_i and 
of the operator (the projection of the total angular momentum J along 
the quantization axis). The associated quantum numbers are ji,j2, ■ ■ ■ ,jn+i 
; J; ki, k2, . . . , kn-i and M, where —J < M < in integer steps. If 

(g) ■ ■ ■ ^W" (g) W"+'^ denotes the factorized Hilbert space, namely the 
(n + l)-fold tensor product of the individual eigenspaces of the (J^)^ 's, the 
operators Kj's represent intermediate angular momenta generated, through 
Clebsch-Gordan series, whenever a pair of J^'s are coupled. As an exam- 
ple, by coupling sequentially the J^'s according to the scheme (■ ■ ■ ((Ji + 

J2) + J3) H \- Jn.+i) = J - which generates (Ji + J2) = Ki, (Ki + J3) = 

K2, and so on - we should get a binary bracketing structure of the type 
(■ ■ ■ {{{W' ®W^)ki ®W'')k2® ■ ■ ■ ® ^^"+')fc„-i)j, where for completeness we 
add an overall bracket labelled by the quantum number of the total angu- 
lar momentum J. Note that, as far as j^'s quantum numbers are involved, 
any value belonging to {0,1/2,1,3/2,...} is allowed, while the ranges of 
the fcj's are suitably constrained by Clebsch-Gordan decompositions {e.g. if 
(Ji + J2) = Ki ^ |ji — 72! < ^1 < 3i +32)- We denote a binary coupled basis 
of (n + 1) angular momenta in the JM-representation and the corresponding 
Hilbert space introduced in (jH)) as 

{ I bi, J2, J3, . . . , Jn+i]^ ; A;2^ . . . , kl_, ;JM),-J<M<J} 

= ni{h) = span{ |b; JM)„} , (4) 

where the string inside [ji, j2, js, . . . , jn+i]^ is not necessarily an ordered 
one, b Bindicates the current binary bracketing structure and the fcj's are 
uniquely associated with the chain of pairwise couplings selected by b. 
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• For a given value of J each Ti^ih) has dimension (2J + 1) over C and 
thus there exists one isomorphism 

ni{h) c^-^+i (5) 

for each admissible binary coupling scheme b of (n + 1) incoming spins. 
The vector space C'^-^^^ is naturally interpreted as the typical fiber attached 
to each vertex v{h) G V of the fiber space structure (0) through the iso- 
morphism (0). In other words, Hilbert spaces corresponding to different 
bracketing schemes, although isomorphic, are not identical since they actu- 
ally correspond to (partially) different complete sets of physical observables, 
namely for instance {J^, J2, J12, J3, J^, Jz} and {J^, Jg, J3, J23, J^, Jz] re- 
spectively (in particular, J^2 and J23 cannot be measured simultaneously). 
On the mathematical side this remark reflects the fact that the tensor prod- 
uct ® is not an associative operation. 

• For what concerns unitary operations acting on the computational 
Hilbert spaces (j3)), we examine first unitary transformations associated with 
recoupling coefficients (Sraj symbols) of SU{2) (j-gates in the present quan- 
tum computing context). As shown in any such coefficient can be split- 
ted into 'elementary' j-gates, namely Racah and phase transforms. A Racah 
transform applied to a basis vector is defined formally a,sTZ : | . . . ( (a 6)^ c) / . . . 
I— *• I . . . (a (6 c)e )/ . . . ; JM) , where Latin letters a,b,c, . . . are used here to de- 
note generic, both incoming {ji 's in the previous notation) and intermediate 
(/cj 's) spin quantum numbers. Its explicit expression reads 

|(a(6c)e)/;M) 

^^(_l)a+5+e+/[(2rf^l)(2e + l)]V2| « ^ ^ I |((a6),c);;M), (6) 

where there appears the 6j symbol of SU{2) and / plays the role of the 
total angular momentum quantum number. Note that, according to the 
Wigner-Eckart theorem, the quantum number M (as well as the angular 
part of wave functions) is not altered by such transformations, and that the 
same happens with 3nj symbols. On the other hand, the effect of a phase 
transform amounts to introducing a suitable phase whenever two spin labels 
are swapped 

\...{ab),...-JM) = {-ir+'-'\...{ba),...;JM). (7) 
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These unitary operations are combinatorially encoded into the edge set E = 
{e} of the graph Gn(y,E): E is just the subset of the Cartesian product 
[V X V) selected by the action of these elementary j-gates. More precisely, 
an (undirected) arc between two vertices v{h) and v{h') 

e(b,b') = (t;(b), t;(b')) e (VxV) (8) 

exists if, and only if, the underlying Hilbert spaces are related to each other 
by an elementary unitary operation (0) or (|7j). Note also that elements in E 
can be considered as mappings (V x C^'^'^^)n — > (V x C^'^+^)„ 

{v{h),Hi{h)) ^ {v{h'),ni{h')) (9) 

connecting each given decorated vertex to one of its nearest vertices and 
thus define a 'transport prescription in the horizontal sections' belonging to 
the total space {V x C^"^+^)n of the fiber space The crucial feature 

that characterizes the graph Gn{V,E) arises from compatibility conditions 
satisfied by 6j symbols in ©, cfr. fn\ . 

The Racah (triangular) identity, the Biedenharn-EUiott (pentagon) identity 
and the orthogonality conditions for 6j symbols ensure indeed that any simple 
path in G„(V, E) with fixed endpoints can be freely deformed into any other, 
providing identical quantum transition amplitudes at the kinematical level. 

To complete the description of the structure (V, C^"^+^, SU{2y)n we 
should call into play M-gates which act on the angular dependence of vec- 
tors in 7i;[(b) by rotating them. The shorthand notation SU{2Y employed 
in (0) actually refers to the group of Wigner rotations, which in turn can 
be interpreted as actions of the automorphism group of the fiber C^"^"*"^. 
Since rotations in the JM representation do not alter the binary bracketing 
structure of vectors in computational Hilbert spaces we might interpret these 
actions as 'transport prescriptions along the fiber'. However, we are going to 
exploit here only computational degrees of freedom associated with j-gates, 
and thus we do not need to explicitate the action of gates along the fiber. 
For this reason, we switch to the notation GniV, E) x C^'^'^^, instead of (j21), 
to represent the computational space of the spin network simulator. 

Let us point out that our model of quantum simulator actually complies 
with a variety of computing schemes, ranging from circuit-type models and 
finite states-automata up to discretized versions of 'topological' quantum 
computation. Inside each of these classes we may also think of different 
encoding schemes to deal with particular problems. In this respect, problems 
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from low dimensional topology, geometry, group theory and graph theory 
(rather than from number theory) turn out to be particularly suitable to be 
addressed in this 'quantum-combinatorial' framework. 

Generally speaking, a computation is a collection of step-by-step transi- 
tion rules (gates), namely a family of 'elementary unitary operations' and we 
assume that it takes one unit of the intrinsic discrete time variable to perform 
anyone of them. In the combinatorial setting described above such prescrip- 
tions amount to select (families of) 'directed paths' in G„(V,i?) x C^'^"^^ 
all starting from the same input state and ending in an admissible output 
state. A single path in this family is associated with a particular algorithm 
supported by the given program. By a directed path V with fixed endpoints 
we mean a (time) ordered sequence 

VL)n = |Vout)n , (10) 

where we use the shorthand notation |vs)„ for computational states and 
s = 0,1,2,..., L{V) is the lexicographical labelling of the states along the 
path. L{V) is the length of the path V and L[V) ■ r = T is the time required 
to perform the process in terms of the discrete time unit r. 

A computation consists in evaluating the expectation value of the unitary 
operator U-p associated with the path V, namely 

(vout I Up I ^in )n- (11) 

By taking advantage of the possibility of decomposing U-p uniquely into an 
ordered sequence of elementary gates, (fTT|) becomes 

(vout I U-P I Vjj^ )„ = [ n (v.+i I Ks,s+i I V, )„, Jp (12) 

s=0 

with L = L{V) for short. The symbol [ Jp denotes the ordered product 
along the path V and each elementary operation is rewritten as Us,s+i {s = 
0, 1, 2, . . . L{V)) to stress its 'one-step' character. 

4 Finite states— automata 

The theory of automata and formal languages addresses in a rigorous way 
the notions of computing machines and computational processes . If A is an 
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alphabeth, made of letters, digits or other symbols, and A* denotes the set 
of all finite sequences of words over A, a language C over A is a subset of A*. 
The empty word is e, the concatenation of two words u and v is simply de- 
noted by uv. In the sixties Noam Chomsky introduced a four level-hierarchy 
describing formal languages according to their internal structure, namely 
regular languages, context-free languages, context-sensitive languages and 
recoursively enumerable languages (see e.g. 11^). The processing of each 
language is inherently related to a particular computing model. Here we are 
interested in finite states-automata, the machines able to accept regular lan- 
guages. A deterministic finite states-automaton (DFA) consists of a finite set 
of states S, an input alphabeth A, a transition function F : S x A ^ S, an 
initial state Sj„ and a set of accepting states Sacc C. S. The automaton starts 
in Sin and reads an input word w from left to right. At the i — th step, if the 
automaton reads the word Wj, it updates its state to s' = F{s,Wi), where s 
is the state of the automaton reading Wi. The word has been accepted if the 
final state reached after reading w is in Sacc- 

In the case of non-deterministic finite states-automaton (NFA) the transi- 
tion function is defined as a map from S x A in P{S), where P{S) is the 
power set of S. After reading a particular symbol the transition can lead 
to different states according to some assigned probability distribution. If an 
NFA has n internal states, for each symbol a & A there is an n x n transition 
matrix for which {Ma)ij = 1 if and only if the transition from the state 
i to the state j is allowed once the symbol a has been read. 
From a general point of view, 'quantum' finite states-automata (QFA) are 
obtained from their classical probabilistic counterparts by moving from the 
notion of (classical) probability associated with transitions to quantum prob- 
ability amplitudes. Computation takes place inside a suitable Hilbert space 
through unitary matrices. Following the measure-once quantum finite- 
automaton (MOQFA) is defined as a 5-tuple M = {Q, S, 6, qo, q/), where: Q 
is a finite set of states; E is a finite input alphabet with an end-marker sym- 
bol (5 : Q X S — > Q is the transition function; 5(q, cr, q') is the probability 
amplitude for the transition from the state q to the state q' upon reading 
symbol a; the state qo is the initial configuration of the system, and q/ is 
the accepted final state. For all states and symbols the function 6 must be 
unitary. At the end of the computational process the automaton measures 
its configuration: if it is in an accepted state then the input is accepted, 
otherwise is rejected. The configuration of the automaton is in general a 
superposition of states in the Hilbert space where the automaton lives. The 
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transition function is represented by a set of unitary matrices U^^cx G S), 
where represents the unitary transition of the automaton reading the 
symbol a. The probabihty amphtude for the automaton M to accept the 
string w is given by 

/Af(w^) = (q/l f/«, |qo) , (13) 

where exphcit form of fKiiw) defines the language C accepted by the au- 
tomaton M. 

The spin network computational space GniV^E) x C^"'^^ (for a fixed n) 
naturally encodes the structures that define (families of) quantum automata. 
Any finite set of binary coupled states belonging to the computational Hilbert 
spaces dH) may represent states of some automaton, while combinations of 
the unitary operations introduced in Q and ((Tj) acting on such states are 
actually transition functions and the amplitude (jl3p complies with expec- 
tation values introduced in ()12|). The inherently step-by-step character of 
transition functions is related to the existence of an intrinsic discrete-time 
variable, denoted by r in the previous section. All such features make the 
spin network the ideal candidate for handling with finite-states (discrete- 
time) automata able to accept any 'quantum language' compatible with the 
algebra of SU{2) angular momenta. In the next section we shall provide an 
explicit example for this class of quantum language. 

5 Quantum automaton calculation 
of link invariants 

As we shall recall below, (colored) link invariants can be obtained as expec- 
tation values of Wilson-type operators in the 3D quantum SU (2) Chern- 
Simons theory for each fixed value of the coupling constant k (related to 
the deformation parameter of the universal enveloping algebra Uq{sl{2)), q = 
root of unity) [Hj. 

From a purely algebraic viewpoint, the basic ingredient for building up such 
link invariants is the 'tensor structure' naturally associated with the repre- 
sentation ring of the Lie algebra of any simple compact group. In the case of 
SU{2) this structure is provided by (tensor products of) Hilbert spaces sup- 
porting irreducible representations together with unitary morphisms between 
them: these objects are collected in the so-called Racah-Wigner tensor cat- 
egory introduced in section 3 [12]. In order to deal with (planar diagrams 
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of) links we need also to specify the eigenvalues of the braiding matrix to be 
associated with the crossings of the links and this is easily achieved by 'braid- 
ing' the Racah-Wigner category. In the present context, it is natural to take 
advantage of quantum group techniques in order to 'split' phase transforms 
like that in ((7j) by assigning different wheights -depending on the deforma- 
tion parameter q- to right and left handed twists. From the combinatorial 
viewpoint, this generalization corresponds to replace the spin network com- 
putational space G„(V, E) x C^"^^^ (encoding the Racah-Wigner category) 
with its g-braided counterpart (G„(V, E') x C^'^'^^) x Z2, where 6j symbols 
in the Racah trasforms © become g-deformed. 

As a general remark, let us point out that the link invariants we are going 
to address are 'universal' in the sense that historically distinct approaches ( 
R-matrix representations obtained with the quantum group method, mon- 
odromy representations of the braid group in 2D conformal field theories, 
the quasi tensor category approach by Drinfeld and the ?>D quantum Chern- 
Simons theory, c/r. ^Hl for a review) are indeed different aspects of the 
same underlying algebraic structure. In this sense the Chern-Simons setting 
describes the universal structure of braid group representations shared by 
all the models quoted above, and in particular the 'universal' (colored) link 
polynomial essentially coincides with the Reshetikhin-Turaev invariant [T7] . 
However, when the group is SU{2), we may speak of 'extended' Jones poly- 
nomials since the Jones polynomial jS] is recovered by selecting the j = \ 
representation on each of the link component (or on each strand of the as- 
sociated braid). On the other hand, the topological quantum field approach 
calls into play links embedded in three-dimensional manifolds, and thus, for 
instance, the invariants can be naturally interpreted as quantum invariants 
of hyperbolic manifolds (complements of links in the 3-sphere). 

Recall that the SU (2) Chern-Simons action for the 3-sphere is 



where A is the connection one-form with value in the Lie algebra of the gauge 
group SU{2) and k is the coupling constant (constrained to be a positive 
integer by the gauge-invariant quantization procedure). The observables of 
the associated quantum field theory are Wilson loop operators defined, for 
an oriented knot K carrying a spin-j representation, as 



kS{A) = ^[ tr{AdA+lAAAAA) 



An Jss 3 



(14) 




(15) 
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and, for a link L made of s components Ki, each labelled by a spin (coloring), 



Wnn..,s [L] = n lU (16) 



1=1 



The expectation values of the above operators are formally given by 
V. \n - m W,,..,^ [L] e-^^(^) 

where the functional integration is over flat S'f/(2)-connections on and 
the generating functional in the denominator will be normalized to 1. Such 
expectation values depend on the isotopy type of the oriented link L and on 
the set of representations labelled by (ji, ...,js) associated with the link com- 
ponents (and do not depend on the metric of the ambient manifold) . To eval- 
uate explicitly such functional averages we have to go through the relation- 
ship between Chern-Simons theory on an oriented, compact three-manifold 
with boundary and the Wess-Zumino-Witten (WZW) conformal field theory 
which lives on the 2-dimensional boundary components. If S*^^\ . . . , S*^^^ 
are the boundary components, (a finite number of) Wilson lines (carrying la- 
bels ji) intersect E*^*-* at some points (punctures) P/*"*. According to the 
axioms for TQFT, we associates with each boundary component S*-*) a finite- 
dimensional Hilbert space denoted by H^^\ In this case the Chern-Simons 
functional (represented geometrically by a cobordism between incoming and 
outgoing boundaries and algebraically by a functor between the associated 
Hilbert spaces) is a state belonging to the tensor product of the boundary 
Hilbert spaces. It can be shown that the conformal blocks of the SU{2)i 
WZW field theory on the boundaries with punctures actually provide basis 
vectors for the Hilbert spaces if*^*-* (the level i of the WZW model is related 
to the deformation parameter q according to g = exp{— 27r2/(£ + 2)}, and in 
turn £ is related to the coupling constant k{> 3) of the CS theory in the bulk 
by £ = - 2). 

Suppose that the oriented hnk appearing in the the expectation value (fTTj) 
is endowed with a plat presentation [TH], namely is actually the closure of 
an oriented braid with 2m strands. If we remove two open three-balls from 
5^ we get two boundaries (both topologically S*^, but with opposite orienta- 
tions) and we can always accomodate 2m 'unbraided' Wilson lines carrying 
labels ji, ^2, ■■■,j2m starting from the incoming (lower) boundary and ending 
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into the outgoing (upper) one. Following the formulation due to Kaul jTH], 
let us start by denoting this 'identity' oriented colored oriented braid as 



where ji = {ji,ei) i = 1,2,... 2m represents the spin ji together with an 
orientation ei,i = ±1 for a strand going into or away from the boundary, 
while stars over the symbols stand for the opposite choice of the orienta- 
tion. A generic colored oriented braid on 2m strand can be generated from 
this identity braid by applying a suitable 'braiding' operator B (written in 
terms of generators Bi, B2, -B2m-i) starting from the lower boundary. The 
resulting colored oriented braid is denoted by 



where labels have been ordered according to the requirement of having a plat 
presentation for the associated link. According to our previous remark, we 
must now select bases in the two boundary Hilbert spaces H^^^ and if^^-* and 
then built up the Chern-Simons functional which will belong to H^^^ ® H^'^\ 
There are basically two possible choices of such bases on the incoming Hilbert 
space H^^^ , corresponding to two sets of correlators in WZW theory involving 
2m primary fields with angular momentum assignments (ji, . . . ,j2m), namely 

I ((jl j2) ■ ■ ■ (j2m-l j2m))[p;r] )] I ((jl (j2j3) • • • (j2m-2j2m-l)) j2m) [t;s] )• (20) 

The shorthand notations [p; r] and [t; s] represent the strings of interme- 
diate spin variables arising from the given pairwise couplings of the in- 
coming (jl, . . . , j2m), {po,Pi,...;ro,ri,...} and {to,ti, . . . ; so, Si, . . .} respec- 
tively. Moreover, the superscript ^ corresponds to the requirement that all 
such pairings must generate SU{2) singlets (spin-0 representations). Uni- 
tary transformations between the two sets of vectors (^01) are implemented 
by the so-called duality matrices, involving g-deformed 6j symbols (c/r. the 
appendix of for their explicit expressions). 

When the outgoing Hilbert space is considered (corresponding to the bound- 
ary S*^^^ endowed with the opposite orientation with respect to S^^^) we have 
to introduce bases which are dual with respect to (j^Uj) . suitably normalized 
according to 

( ((jl J2) • • • (j2m-l J2m))[j,;r] I ((jl J2) • • • (j2m-l J2m) ) [|,';r'] ) = ^(p)(p') ^ir){r') 




(18) 




(19) 



(21) 
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and 



(((Jl(j2j3)(j2m-2j2m-l))j2m)°;s] | ((jl (7273) • • • ( j2m-2 j2m-l ) ) j2m) ° ';^'] ) = 5(t)(t') S{s){s')- 

(22) 

The basis vectors (j^n|) are eigenf unctions of the odd and even braid genera- 
tors, respectively 

52i+l|((jlj2)---(j2«+lj2«+2)--.)b;r]) = ^p, (j2«+l , j2i+2) | ( ( jlj2) • • • (j2i+2j2i+l ) • • •) fpir] ) 

B2I I ( (j2m jl ) • • • (j2« j2/+l ) ■■■)l;s]) = K U2I , J2«+1 ) | ( ( j2mjl ) • • • (j2«+l J2« )--)[t-s]) 

(23) 

with eigenvalues given by the expressions 

A, (j, j') = AW(j>j') = (-)^'+^"-^g('=^+'=^')/2+c„,i„(,,,,)-c./2 ^ 

K 0.3') = i>^i'\j,f)r' = if ee' = -1 (24) 

where q is the deformation parameter, z G {pz,^;} and = z{z + 1) is the 
quadratic Casimir for the spin z representation. Thus A^"*"-* ( j, j') is the eigen- 
value of the matrix which performs a right handed half-twist in strands with 
the same orientation, while is the eigenvalue of the matrix which 

performs a right handed half-twist in strands with the opposite orientations. 
The expectation value of the Wilson operator ()17p for any link L presented 
as the plat closure of an oriented braid over 2m strands can be recasted into 
the expectation value of the braiding operator B associated with the oriented 
braid (fT^ according to 



Vj^jl-M W = n [2J. + 1] 



X 



^i=l 



((('.':)■■■ («.));„,! «(f; I : t tj|((-o ■■■(-&));„,) 

(25) 

This expression (which is normalized according to the standard conventions) 
gives the value of the extended Jones polynomial of the link L. The oper- 
ator B is expressed in terms of (a finite sequence of) braid generators ()23|) . 
suitably converted in terms of the current basis vectors by acting with g- 
Racah transforms on the even subset, while [2jj + 1] is the g-dimension of 
the representation jj. 
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Coming back to the definition of tlie quantum finite state-automaton 
(QFA) given in the previous section, namely a 5-tuple M = {Q, S, 6, qo, F), 
we recognize that, given a plat presentation of a link L, we can always build 
up a QFA Ml which recognizes the language of the braid group and whose 
associated probability amplitude (fT^ is given by the value of the extended 
Jones polynomial (j25|) . More precisely: qo is a particular ket vector of the 
type (j2ni); S is the alphabeth given by the (2m — 1) braiding operators (j2Sl); 
F is a set of final stated suitably selected among the bra vectors in (pUj). 
Note that, according to the expression given in ()25|). F must contain (sin- 
glet) states which possess the same pairwise coupling structure as the initial 
state qo exhibits. Owing to the topological properties of the plat presenta- 
tion, this means that final states may differ from qo by a permutation on the 
string (jij2 • • -jm) and thus we can actually build up a family of m! automata 
out of one initial state qo- Once a final state q/ has been selected (the right 
permutation can be singled out in a fast way even by a classical machine) 
the evaluation of is carried out by the automaton linearly in the length 
of the 'word' B. The length of B is the number of { braiding operators -|- 
g-Racah transforms } entering into the explicit expression of B and equals 
the number of crossings of the plat presentation of L. It can be easily rec- 
ognized that such an automaton can be associated with a particular path in 
the g-braided spin network computational space (G„(y, -E) x C^"^+^) x Z2 
for n = 2m — 1 since the process carried out here is a particular realization 
of the general computing procedure described at the end of section 3 (c/r. 

What we have really done is to 'encode' the combinatorial structure un- 
derlying quantum SU{2) Chern-Simons theory (and the associated WZW 
boundary theory) at some fixed level i into the above abstract g-braided 
S'[/(2)-decorated graph. This does not mean, of course, that we have got a 
quantum algorithm in the proper sense, since the encoding map could not be 
'efficient' (nor efficiently approximated) with respect to standard models of 
computation (Boolean circuit, Turing machine). We are currently addressing 
such open problems. 
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